Wave packet Rabi oscillations from phase-space flow in mesoscopic cavity QED 
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The dynamics of the Jaynes-Cummings Hamiltonian is studied in the semiclassical limit by re- 
ducing it to quantum normal form and deriving the phase space dynamics of classically smooth 
observables. The resulting systematic semiclassical approximation is used to analyze the Rabi os- 
cillations of a photon wave packet revealing subtle consequences of the uncertainty principle. 
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Introduction The Jaynes-Cummings (JC) Hamilto- 
nian of a two-state system coupled to a harmonic os- 
cillator is a basic model of quantum optics and the the- 
ory of atom-light interactions. This simple Hamiltonian 
describes quite accurately the quantum electrodynamics 
of strongly coupled Rydberg atoms in ultra-high Q mi- 
crowave cavities [1]. Because of the universal nature of 
the JC Hamiltonian, it models many other quantum sys- 
tem with interactions between discrete and continuous 
degrees of freedom, like hyperfine levels in trapped ions 
2] and spin-orbit interactions in semiconductors [3| . It is 
also a building block for more complicated systems like 
micromasers [J] and electromagnetically induced trans- 
parency cavities [a]. 

In the rotating wave approximation (RWA), where 
non-resonant terms in the Hamiltonian are neglected, the 
JC model becomes integrable, with a conserved quan- 
tity in addition to energy — the polariton number that 
gives the total number of excitations, field and atom. Po- 
lariton number conservation makes it easy to diagonalize 
the Hamiltonian, since the polariton eigenspaces are only 
two-dimensional, and makes the dynamics of states with 
a small photon number simple. Nevertheless, the dynam- 
ics of states that have many significant energy-state com- 
ponents, including the important example of phase-space 
localized wave packets, can show considerable complexity 
[a], including squeezing pj, wave packet Rabi oscillation, 
collapse-revival [8|], and entanglement [3]. 

Recently the interest in wave packet dynamics in the 
JC system has increased, as it has become a subject of 
direct experimental investigation in cavities excited to 
mesoscopic coherent states with a moderate number of 
photons [lfj • The experimental systems are of course also 
subject to dissipation and decoherence, but the coherence 
lifetimes that have been achieved ensure that the actual 
dynamics can be approximated by that of an ideal system 
during the initial time interval of wave packet evolution, 
where the complex dynamical phenomena occur. 

Complex wave packet dynamics can be exhibited by 
the integrable JC system because its energy spectrum is 
a nonlinear function of the number of polaritons. Many 
of the complex dynamical phenomena occur also in the 
simpler case of nonlinear oscillators 111, [l2| , where they 
can be analyzed using semiclassical nonlinear dynamics. 
Here we show that an analogous semiclassical analysis is 



applicable also to the JC dynamics. 

As always, the semiclassical approximation is valid 
when the typical action scales are significantly larger than 
Planck's constant. In the JC system this condition holds 
either when the intrinsic action scale Br (denned below), 
that is a measure of the detuning, is much larger than h, 
or when the initial state is sufficiently excited energeti- 
cally. Then the field degree of freedom becomes semiclas- 
sical but, unlike nonlinear oscillators, the JC system also 
has the atomic polarization internal degree of freedom 
that remains fully quantum mechanical. This contrast 
is reflected by the much faster dynamics of the polariza- 
tion versus the field: The ratio of the natural frequencies 
of the polarization and the field dynamics diverges in 
the semiclassical limit. Our semiclassical analysis, that 
retains an exact quantum-mechanical treatment of the 
internal degree of freedom, is in this way different from 
that of cavity dynamics of an ensemble of atoms [13| . 

Wave packet propagation studies [a l9( have provided 
a semiclassical analysis of the JC system maintaining a 
quantum internal degree of freedom, but their utilization 
of the classical dynamics is limited because Scrbdinger 
picture wave packet propagation does not have a direct 
classical limit. Here, in contrast, the dynamics is ob- 
tained by evolving classically smooth observables in the 
Heisenberg picture. The dynamics of the observables is 
cast in terms of their semiclassical phase-space quasi- flow 
111 . ll4J , defined as their time-dependent Weyl represen- 
tation. As a consequence of the presence of internal 
degrees of freedom the quasi-flow of the JC system is 
matrix-valued, giving the components of an operator in 
polarization subspace at each phase-space point. 

The focus of this work is the Rabi oscillations of wave 
packets. The semiclassical quasi-flow leads to the explicit 
calculation of the Rabi dynamics beyond the leading- 
order classical substitution, demonstrating subtle quan- 
tum effects in the interaction between the slow and fast 
degrees of freedom: The entanglement of eigenstates re- 
sults in the polarization oscillations shown in Eq. ([6]) and 
Fig. [TJ while the quantum field uncertainty causes the 
field Rabi oscillations shown in Eq. {7J and Fig. [2] 

The polariton quasi-flow The JC model describes a 
single electromagnetic mode of frequency u>, modeled by 
a harmonic oscillator with a lowering operator a, inter- 
acting with a two-state system of energy gap hv with 



a lowering operator a. Under the RWA, where the non- 
resonant terms aa and a' a' are neglected, the JC Hamil- 
tonian is 



f(B + H)b that holds for an arbitrary function 



H = huj 



a) a + aa^ 



hv- 



■ + %(a t cr + acr t ) . (1) 
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When the two-level system models an atom, g is an ef- 
fective dipole interaction. The commuting a and a op- 
erators obey Bose and Fermi algebra rules, respectively, 
i.e., [a,aJ] — 1, a 2 — 0, [a, a^] + = 1. 

The RWA JC Hamiltonian is integrable because it 
commutes with the polariton number operator n po \ = 
a^a — a^a. In other words, the system excitations are po- 
laritons 15| with an internal two-state polarization de- 
gree of freedom. For the purpose of studying dynamics 
we therefore define the (single-mode) bosonic polariton 

annihilation operator b = ^fh{a^ a + <j<j^ J f^r) a ' an( i the 
fermionic polarization lowering operator r = era 1 } 



that commutes with it. That is, b< and b create and 
annihilate (respectively) polariton quanta, and r' and r 
raise and lower (respectively) the polarization. We define 
also standard SU(2) generators for the internal degree of 
freedom of the polaritons, n = r 4 

rt 



r 2 = i(r - r 1 *) 

t- 3 = T I T — TT \ 

In terms of b and r the Hamiltonian takes the form 
H = ujB + St 3 + XvBtx , where 5 = | (u — w), A = y/hg , 

B = btf = h(n po \ + 1), and V B is the unique positive 
operator whose square is B. The advantage of the po- 
lariton representation is that the Hamiltonian depends 
on the polariton field only through the combination B; 
that is, the Hamiltonian has been reduced to its quan- 
tum normal form [16([l8|. At this stage B can be han- 
dled as if it were a c-number, letting the Hamiltonian 
be further reduced by invariant polarization space de- 
composition into H = P+(B)e+(B) + P-(B)i-(B), in- 
troducing the definitions s±(B) — iqB ± e(B), e(B) = 
Av B + B R , where B R — {S/X) 2 is the natural action 
scale of the system. e±(B) are scalars in the internal 
degree of freedom space, i.e. they commute with r and 
rt. P±(B) = |(1 ± (s(-B)ti + c(B)t 3 )) are orthogonal 
projections such that P + + P_ = 1, P+P- = P-P+ = 0, 



with c(B) = JB R /(B + B R ), and s = Vl-c 2 . The 

P±(B) are projections on invariant dressed polarization 
subspaces that, unlike the bare polarization projections 
i(l ± t 3 ), do not commute with the field b. 

The reduction of the Hamiltonian to its normal form 
simplifies the solution of the Heisenberg equations of mo- 
tion. The dynamics of the polarization operator T3 is 
given by 

r 3 ,t = (s(B)n + c(B)e T3 )c(B) - (3(3)0(3)^ - r 3 ) 

x cos(Q(B)t) + s(B) sm(n(B)t) , (2) 

where fi = %e is the Rabi frequency. The dynamics of 
the polariton field b is derived using the identity bf(B) = 



b t = (P+(B) + P_(B)e- l0 ^*)P + (B + h)i 



-iun{B)t 



(3) 



+ (P + (B)e*°^* + P_(P))P_(P + hy^^e-^b , 

defining the (field-dependent) phase frequency ui% — 
i(e(- + K) — £(•))• &h is semiclassically slow and has a 
well-defined value w = e' in the classical limit h — > 0. 

The dynamical analysis presented so far is exact, but 
the conversion of Eqs. (O and © to physically measur- 
able quantities requires taking matrix elements with re- 
spect to the initial state, and for this purpose the semi- 
classical approximation will be assumed, letting H — ► 
keeping other system parameters fixed. The analysis will 
rely on the classical smoothness of T3,t and bt, and this 
requires the time t to be of o(l) in the semiclassical limit. 
As shown below, the time of validity encompasses many 
Rabi oscillations and their collapse but not their revival. 

The semiclassical approximation will be made in the 
Wigner-Weyl representation of the JC system that is a 
special case of the standard Wigner-Weyl representation 
of many-component systems 17 1 , in itself a straightfor- 
ward generalization of the scalar Wigner-Weyl represen- 
tation. Because the polarization degree of freedom re- 
mains quantum-mechanical, the symbol range consists of 
the operators in the two-dimensional polarization Hilbert 
space. In the same manner, the Wigner "function" W is 
also operator- valued and the expectation value of an ob- 
servable is 



L 



d 2 a 

~2T 



tr L(a, a)W (a, a) 



(4) 



where L is the Weyl representation of the operator L, a 
is a complex phase-space coordinate, and tr stands for 
trace on the two states of the internal degree of freedom. 
The dynamical expectation values presented below 
were calculated using Eq. Q after deriving the quasi- 
flow, i.e. the time dependent Weyl representation, of 03 
(= t 3 ) and of a from Eqs. ^ and (J3J). Because these 
operators are classically smooth, the leading term quasi- 
flow is obtained simply by substituting the classical field 
amplitude a and its complex conjugate (respectively) for 
the operators b and tf (note that the difference between 
the a and b operators is of higher order in K) . The higher 
order terms were obtained using the standard rules of 
the Moyal calculus [T7[ . The straightforward but techni- 



cal details of this calculation are omitted. 

Semiclassical Rabi oscillations We proceed to a con- 
crete demonstration of the results by the Rabi oscilla- 
tions of phase-space localized pure initial states. Specif- 
ically, the initial states are taken to be separable in the 
bare field-atom representation with a coherent field state. 
This means that the initial Wigner "function" is of the 
form 2wexp(— 2\a — cto| 2 ), where w is a projection oper- 
ator in the atomic two-dimensional Hilbert space. The 



Wigner function is combined with the quasi-fiows to cal- 
culate physical expectation values, by carrying out the 
Gaussian phase-space integral and the two-dimensional 
trace operation in Eq. (fj|. 

Consider first the case where the atom is prepared in 
the excited state, w = 5(1+0-3), and the observable to be 
measured is the mean polarization (03) that gives the dif- 
ference between the probabilities that the atom is in the 
excited state and in the ground state. In a classical field 
of intensity equal to the mean intensity of the quantum 
field, the two-level system would undergo periodic Rabi 
oscillations of frequency Q(Aq), where A = h(\ao\ 2 + 5) 
is the expectation value of the action in the coherent ini- 
tial state |ao). However, because of the field uncertainty, 
the Rabi oscillations have a finite lifetime [8j after which 
they collapse. The semiclassical calculation agrees with 
this result: 

(<7 3 t) = c(Ao) 2 + s(Ao) 2 co S (n(A )t)e-^ A ^ n ^ A ^ 2 

(5) 
exhibits Rabi oscillations of frequency fl(Ao), super- 
exponentially damped with the rate \ZKAqH' (Aq) . The 
collapse time is therefore asymptotic to ftr 1 ' 2 ; in particu- 
lar it diverges in the classical limit as expected. It should 
be mentioned that Eq. ([5]) is obtained from Eq. (F2J) us- 
ing only the leading classical substitution term in the 
quasi-flow; therefore, the same collapse dynamics would 
be observed with a random classical electromagnetic field 
with a Gaussian distribution. 

A more subtle consequence of the quantum uncertainty 
in the electromagnetic field is the impossibility of con- 
centrating a state that is separable in the bare atom- 
field decomposition on one of the ± dressed polarization 
subspaces, because the energy eigenstates are entangled. 
As a result, the dynamics of any separable initial state 
exhibits Rabi oscillations between the two internal po- 
lariton states. The separable state will be maximally 
localized on the + subspace if we choose w = P + (A ). 
The mean polarization obtained with this initial state is 

(a 3t ) = c(A ) - n^^s(A )(A 9'(A Q ) (6) 

+ is(A )) 8 in(n(A))t) exp(-iM (fi'(A )i) 2 ). 

Beside the constant term, that is the only one to survive 
in the classical limit, the polarization displays Rabi oscil- 
lations with an amplitude proportional to H/Aq, that is 
growing in time. The amplitude does not grow without 
bound, however, because it is also subject to collapse. 
Since the collapse time grows as hr 1 ! 2 in the semiclassi- 
cal limit, the peak amplitude is of Oifi 1 ^ 2 ). A comparison 
of the semiclassical calculation with results obtained by 
direct numerical solution of the Schrodinger equation for 
a particular parameter choice is shown in Fig. [TJ showing 
good agreement for a moderate number of photons. 

A different manifestation of the quantization of the 
electromagnetic field is the oscillations of the field ampli- 



tude (a) with the fast Rabi frequency. One cause for these 
oscillations is the periodic energy exchange between the 
field and atom. More surprising, though, are the Rabi os- 
cillations in the polariton field, caused by the non-trivial 
commutator [b, b'] that generates the off-diagonal terms 
in Eq. ©. 

As an example consider again the case where the atom 
is initially in the excited state, w = ^(1 + t 3 ). In the 
frame rotating with the optical frequency u), the expec- 
tation value of the field amplitude is 

(a t ) e lu}t = (cos(w(A)*)) - ic(A ) sin(w(A,)t)) (7) 

-^ S (Ao) 2 (l-e-^°)*e X p(-iMo(n't) 2 ))ao. 

The terms in the second line express the Rabi oscillations 
of amplitude proportional to h/ Aq. The terms on the first 
line, on the other hand, result from weighted average of 
the two diagonal elements in Eq. ([3]) that describe adi- 
abatic phase oscillations in the ± dressed polarization 
invariant subspaces frequency ±u) (respectively). In the 
classical substitution approximation these phase oscilla- 
tions lead to the splitting of the diagonal components of 
the Wigner function [g, Q, LLC} • 
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FIG. 1: (Color online) The mean polarization of an initially 
coherent photon state, and atomic polarization optimally lo- 
calized on the '+' invariant subspace, with Br = 6.25ft. The 
continuous blue line is derived from the semiclassical approx- 
imation Eq. (JS| and the black dots from a direct numerical 
solution. Good agreement is observed already in the top panel 
where the mean polariton number is 8 while the bottom panel 
with 50 polaritons shows agreement to less than 3%. The rel- 
ative error term in the semiclassical approximation is 0(H 1 ' 2 ). 

The time evolution of the the field amplitude as ob- 
tained from Eq. ([7]) is shown in Fig. [5] for a particular 
parameter choice, and compared with the values obtained 
from the numerical solution. The "baseline" of the Rabi 
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FIG. 2: (Color online) The expectation of the real quadrature 
of the field amplitude of an initially coherent photon and ex- 
cited atomic state with Br = 6.25ft. The continuous blue line 
and the black dots have analogous meanings to those of Fig. 
[U while the dashed red line shows only the adiabatic terms 
in the semiclassical field dynamics Eq. (0. Good agreement 
is observed for the lower polariton numbers, 4 (top) and 8 
(bottom) since the relative semiclassical error is 0(h). 



This result is the starting point for the analysis of quan- 
tum electromagnetic effects as an expansion in powers of 
h. The concept was demonstrated using photon coher- 
ent states, but clearly holds for any phase-space local- 
ized initial state. Furthermore, the idea is not limited 
to the analysis of Rabi oscillations: It is applicable to 
most other quantities of interest in the JC system such 
as squeezing and entanglement, as well as to other multi- 
component quantum system including some cited in the 
introduction. 

The approximation of quantum electrodynamic phe- 
nomena by classical ones is, like all semiclassical dynam- 
ical approximations, nonuniform in time. In our case, 
the time scale of validity is set by the Heisenberg time 
tfi = (ft^') _1 j where the discreteness of energy levels be- 
comes manifest with the revival of oscillations 8], that 
is not perturbative in the quasi-flow. Neither is the 
semiclassical quasi-flow well-suited to study quantum in- 
terference phenomena. These phenomena are of course 
amenable to semiclassical analysis, but they require the 
use of more elaborate tools and considerably heavier anal- 
ysis, in contrast with the purpose of this work. 



oscillations is provided by the adiabatic terms, whose 
contribution is shown in the figure with a dashed line. 
The adiabatic motion is so slow, that in the time interval 
before the collapse of the Rabi oscillations, only a frac- 
tion of a phase oscillation is achieved. However, because 
of the much larger amplitude of the adiabatic motion, 
its effect is comparable to that of the fast Rabi oscilla- 
tions, illustrating the singular nature of the semiclassical 
limit, where phenomena occurring on widely separated 
time scales have to be considered simultaneously. 

Conclusions The dynamics of quantum systems, in- 
cluding integrable ones, becomes computationally com- 
plex for phase-space localized wave packets, or other 
initial states having many nonzero energy components. 
However, this difficulty can be bypassed when expecta- 
tion values of classically smooth observables are sought, 
and the utility of this approach for multi-component sys- 
tems was demonstrated with the semiclassical analysis of 
the Rabi oscillations in the Jaynes-Cummings system, a 
phenomenon that is a direct manifestation of its quantum 
multi-component nature. 

Computability is not a practical concern for the JC sys- 
tem, but the semiclassical analysis provides a thorough 
and systematic theory of the dynamics that is lacking in 
pure numerical analysis. The results presented in Eqs. 
©-([7]) illustrate the fact that the Rabi oscillations of a 
phase-space localized wave packet are approximated in 
the leading term by oscillations in a classical field with 
parameters determined by the center of the wave packet. 
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